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Using a microscopic model for stochastic transport through a single quantum dot that is modified 
by the Coulomb interaction of environmental (weakly coupled) quantum dots, we derive generic 
properties of the full counting statistics for multi-stable Markovian transport systems. We study 
the temporal crossover from multi-modal to broad uni-modal distributions depending on the initial 
mixture, the long-term asymptotics and the divergence of the cumulants in the limit of a large 
number of transport branches. Our findings demonstrate that the counting statistics of a single 
resonant level may be used to probe background charge configurations. 
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INTRODUCTION 



II. ILLUSTRATIVE PICTURE 



The coexistence of several stationary states for a given 
set of parameters is typically referred to as the phe- 
nomenon of multi-stability. Multi-stable behavior is 
found in a wide variety of systems in different disciplines 
of science, as e.g. biology^, chemistry^, neuroscience^, 
laser physics^, and semiconductor physics^. 

In transport systems, multi-stability is characterized 
by the existence of more than two distinct branches in 
the transport characteristics with hysteresis and switch- 
ing in between. Some prototype examples for correspond- 
ing electronic systems are superlattices^, double-barrier 
resonant tunneling diodes^, and nano-electromechanical 
systems^. If the transport is entirely governed by 
stochasticity, e.g. as in single-electron transport^, the 
current alone might not reveal the multi-stable character 
and other more sensitive tools are required. As has been 
shown for bistable systems^S, the counting statistics^i 
may serve as such a tool. Recently, in Ref.— the measure- 
ment of a bi-modal distribution of quantum dot tunnel- 
ing has indicated the interplay of fast and slow transport 
channels not visible in the current. 

In this work, we present a generic approach to study 
Markovian transport systems with multi-stable behav- 
ior. Starting from a microscopic model for one trans- 
port channel with an environmental control system we 
derive a Master equation for Counting Statistics with an 
arbitrary number of transport branches. The resulting 
Liouville super-operator has a simple and scalable block- 
tridiagonal structure. Even though there exists a unique 
steady state, the counting statistics and higher-order cu- 
mulants display clear signatures of multi-stability such 
as multi-modal or broad distributions and diverging cu- 
mulants. We provide results for the temporal evolution 
and long-term asymptotics of the statistics and discuss 
the limit of a large number of coexisting current branches 
analytically. We emphasize that our approach is not re- 
stricted to electrons as transferred entities - in principle, 
stochastic multi-stable transport of any countable object 
can be studied by this means. 



A single transport channel (single resonant level, quan- 
tum dot) influenced by k background charges distributed 
on a collection of N sites will experience an effective shift 
of its charged energy state (e^^ + kU , where U is 

the Coulomb repulsion). Attaching two reservoirs {L, R) 
held at different chemical potentials - see Fig. [Ija) - will 
now induce transport through the dot at rate F, which 
strongly depends on the number of background charges. 
Effectively, this leads to shifted currents in the transport 
channel 



h. = - 



fdea + kU)- fR{ed + kU) 



(1) 



where fL,R{oj) = [l + e+^^'^'F^/^)] denote the Fermi 
functions of the respective reservoir with bias voltage V , 
see Fig.[TJb). Coupling the background charges to differ- 
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FIG. 1: (Color online) (a) Scheme of the model, (b) Re- 
gion of multi-stability: Current-voltage (Transport) charac- 
teristics for N — 7 control dots; thin red curves: A'' -f 1 indi- 
vidual currents Ik (fc = 0, 1, . . . , 7), thick blue line: average 
current (I) which is actually observed; Borders of colored re- 
gions: standard deviation {{P)) /2 of the current which 
becomes divergent for 77 — > in the multi-stability region. 
Inset: Broad uni-modal Distribution function log(P)/t vs. 
current for large times and V* =7 (thick blue curve); distribu- 
tions for individual currents 1^ (thin red curves). Parameters: 
N =7, ed=3, /3=3, f/=0.15, 7^=0.001. 
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ent reservoirs L, R with rate ry ^ F will cause slow ran- 
dom switching between the different transport channels 
Ik (for two currents fc = 0, 1 known as random telegraph 
noised). When the lead temperature is comparable to 
the Coulomb interaction « J7, this gives rise to a pro- 
nounced region of multi-stability around V* « 2ed + NU 
in the current-voltage characteristics, see Fig. [IJb). In 
this region, the second-order cumulant of the current 
{{l^}) diverges for 77 as indicated. The inset shows 
the corresponding broad long-term distribution of cur- 
rents at V* in comparison with the distributions of indi- 
vidual currents. 



III. MICROSCOPIC MODEL 
A. Hamiltonian 

We consider the total Hamiltonian 

N N N 

i—l j—i+l 



i=l 



N 



Uidld^d^'d + ^ iOkacl^Cka 
i—l ka 

■^tka dcl^+CkadA 



N 



d^cl^ + Ckad\ 



(2) 



ka i—l 



where d, di, and Cka annihilate electrons on the trans- 
port dot, the i"^ control dot, and the mode k on lead a 
(with energy lUka), respectively. In addition, we consider 
the symmetrized wide-band limit, where the transport 
dot tunneling rates F = 27r^^ |tfca|^(5(w — Uka) and the 

I -12 

control dot tunneling rates 77 = 27r^^ yka] '^('^ ~ '-^ka) 
become independent of energy and lead. The parame- 
ters Ui denote the Coulomb interaction between trans- 
port and control sites, whereas Uij represents repulsion 
between electrons within the control system. We assume 
that the spectrum of the system Hamiltonian is only near 
but not exactly degenerate Ui ~ U , Uij w Uc, and Ci ~ 0. 
These simplifications are not crucial for the occurrence of 
multi-stability, but rather allow for an analytic treatment 
in the following. 



B. Liouvillian 

We perform our analysis within the Born- 
Markov-secular approximation scheme which can 
be alternativel y-'^^' derived with a coarse graining 
method in the limit of infinitely large coarse graining 
times r. Provided the system energy spectrum is 
non-degenerate and the time scales are larger than 
the inverse minimum level splitting, the Liouvillian 



couples only the diagonals of the density matrix in 
the system energy eigenbasis with each other (see 
alsoiS). Since we are interested in observable effects 
of multi-stability in the current through the transport 
dot, we introduce a virtual detector in the right lead 
via the replacement in the tunneling Hamiltonian 
^"-Ifl ~^ ^^kR ® ^^'^ CkRd^ CkRd^ (g) b, where the 
detector operator = J2n I" + 1) i*^! increases the 
detector outcome each time an electron is created in 
the right transport lead. Treating the tensor product 
of dot and detector Hilbert spaces as the system, we 
arrive at an n-resolved master equation of the form 



{n\p\n) 



Lop' 



(») 



L+P' 



(n-l) 



L.p 



(»i+i) 



which couples different realizations of the dot den- 
sity matrix ~ each valid for different particle num- 
bers n measured in the detector - with each other. 
This coupled system can be further reduced by 
Fourier-transformation p(x,i) = X)n 
where x is the counting field, which leads to 
pix,t) = {Lo + e+'XL+ + e-'XL_) p(x,t) = i(x)p(x,i)- 
Due to the pcrmutational symmetry, it is convenient to 
denote the corresponding eigenstates for TV control sites 

N 



by \N,k,e;nd), where £ £ U,2,. 



k 



arbitrarily 



labels all the configurations with < k < N electrons 
distributed on the N control sites, and € {0,1} de- 
notes the occupation of the transport dot. When we trace 
out the configuration of the control dots for a given total 
number of control charges k by defining the 2x2 matrix 
Pkix^t) = i^'f'^^l Pix,'t)\N,k,£), the LiouviUian in 
this basis assumes for < fc < the form 



Pkix,t) = TC''°\x,ed + kU)pk{x,t) 

+T^{N -k + l)Jtpk-i{x,t) 
+V{k + l)jr'pk+i{x,t), 



r,ClT''°'pk{x,t) 



(3) 



where the newly introduced super-operators arc 2x2 ma- 
trices, which obey J7q™ = and J7^"* = at the bound- 
aries. This defines a 2(A^ + 1) x 2(iV + l)-dimensional 
Liouvillian super-operator L{x) with a block-tridiagonal 
structure. The detailed structure of the reduced Liouvil- 
lian super-operators follows from a rigorous microscopic 
derivation, it may however also be understood from sim- 
ple phenomenological reasoning: 

(i) The multi-stable (fast) part has block-diagonal 
structure, where the N + 1 block matrices correspond 
to the Liouvillian of a single resonant level - shifted by 
the Coulomb interaction with k control charges 



£"°*(x,^) 



e-'Vfl(^) -/i^(^) 



(4) 



where fa{^) = [1 — /a('*^)]- Evidently, when 77 = 0, these 
matrices give rise to the multi-stable currents in Eq. ([T]) . 
Since we have traced out the different control dot con- 
figurations, it also becomes obvious that the associated 
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currents Ik are actually 



-fold degenerate. These 



degeneracies may be lifted (and thereby become observ- 
able) when the assumed symmetries of the Hamiltonian 
are absent. 

(ii) The remaining part of the Liouvillian (which ap- 
pears as slow when rj <i; F) consists of the control system 
jump super-operators 



J out \ 
k = 2^ 

a&{L,R} 



fa{{k-l)Uc) 

fa {{k - l)Uc + U) J ' 



/- {kU, + U) 



(5) 



which depend on the control system occupation number, 
and a trace-conserving part /3g™'™^ = —kj^^\ — {N - 
k)J^i- The scalar coefficients in Eq. ^ arise, since for 
any control configuration with fc — 1 charges there are 
iV — fc -I- 1 different possibilities to obtain a control con- 
figuration with k charges. Similarly, for a configuration 
with k + 1 charges, each single charge leaving the control 
sector constitutes an equivalent jump channeUi. In addi- 
tion, we note that the control jump super-operators ([5]) 
must assume diagonal form in the sequential tunneling 
regime for the basis chosen. 



C. Transport observables 

The probability for obtaining n tunneled particles af- 
ter time t is given by Pnit) = Tr{p(")(t)}. It follows 
that the moments of P„ it) may be directly obtained from 
the Fourier-transformed Liouvillian by suitable differen- 
tiation of the moment generating function (MGF) 



M(x,i)=Tr{e^W*p} 



(6) 



(where Tr {(po, • ■ • , Pw)} = Y^k^aPk)) with respect to 
the counting field x- The initial density matrix p is typ- 
ically chosen as the steady state L{0)p = 0, since one is 
usually interested in long-term cumulants. The matrix 
exponential is significantly harder to evaluate than the 
matrix inverse, such that we consider the Laplace trans- 
form 



Mix,z) = Tr 



1 



zl-Lix) 



(7) 



of the MGF instead. 

For example, the moments of Pn{i) are obtained via 



{n'^){t) = {-id^fM{x.t). 



(8) 



The full distribution, however, is obtained by inverse 
Fourier transform 



P„(t) = ^ / il/(x,i)e-™^dx. 



(9) 
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FIG. 2: (Color online) Current-voltage characteristics for a 
nontrivial non-thermal stationary state (jlOp . (a) for various 
temperatures and A'^ =10; at small temperatures (;9A ^1) 
the number of control dots A'^ can be probed in a current mea- 
surement, whereas at /3A <1 one has to make use of count- 
ing statistics at f < 77^^ (see Fig. [3]). (b) for various A'' and 
/JA =40; for Af -H- 00, f/ ->0, < A and /3A >1 the current 
becomes linear I ~ since transport takes place through 



a continuum with finite band-width A 
the inset). 



IV. RESULTS 



NU (as sketched in 



A. Analytical steady state 

When transport and control dots are coupled to leads 
with the same chemical potential {fa = fa) and = 
as well a.s U — Uc, the steady state of Eq. at x = is 
Cjy (poi ■ ■ ■ Pn) , where the partial vectors read 



P 



Pk = 



k~l ■ 

.£=0 . 



N 



n 2 

l=k+\ 



■Pe 



2-pk 
Pk 



(10) 



with Pk = fL{kU) + fR{kU) and Cn follows from nor- 
malization. The corresponding current / = Tr {L'(0)p} is 
the weighted sum of partial currents Ik ^ ■ The current- 
voltage characteristic exhibits 2N steps for small temper- 
atures /?A ^1 (A = NU) which become smeared out for 
/3A <1 as shown in Fig. [2ja). Therefore, for sufficiently 
low temperatures we are able to probe the number N just 
by a current measurement. At larger temperatures, how- 
ever, this fails and the counting statistics will provide a 
proper tool for that purpose (see Fig. [3] and correspond- 
ing discussion). In the thermodynamic limit {N — > 00, 
U -^0) such that the band-width A stays finite [spectrum 
becomes continuous as sketched in inset of Fig. [2Kb)] the 
characteristic is linear for /3A ^1 and \V\ < A with 
Ohm's resistance of 2A/ eF [compare Fig. El^b)]- 



B. Full Counting Statistics 

The model ([2]) shows very rich behavior. We, however, 
choose some limiting cases to illustrate the multistable 
properties in the following: 
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(i) In the infinite bias limit /l — >■ 1 and /r ^ 0, the 
MGF coincides with that of a single resonant level, for 
which we obtain 



50 100 150 200 250 



Pn>0{t) = 



-rt 



X 



2(2n + 1)! 
X [2n{2n + 1) + (4n + 2)R + (Tt)^] 



(11) 



and Pn<oit) = 0, such that the counting statistics will 
not reveal any multi-stable properties [compare Fig. (Hb) 
for large V]. 

(ii) When the control leads are at infinite bias, i.e. 
fliij) 1 and fjfiuj) such that J^" = J"^""* = 1, 
and the transport leads are at high bias (/^(w) ?» 0), 
one may for sufficiently low temperatures have fLi^d + 
(k < k)U) = 1 and /^(ed + (fc > k)U) = for some 
k 6 {1,...,A^}, which leads to only two different currents 
(bistable case). The detailed form of the Liouvillian and 
its counting statistics then depends on N and k, but the 
whole class of bistable models is amenable to analytic in- 
vestigations. In the simplest case of fc = = 1, we have 
for 77 = a bimodal distribution: Half of the distribution 
follows the evolution of a single resonant level dTTI) . and 
the other half remains localized at 71 = for all times. 
The situation becomes non-trivial for finite ry, which is rc- 
fiected in the recursive relation for the Laplace transform 
P„+i(z) = T{z)Pn{z) for 71 > 2, where F{z) has four dif- 
ferent first order poles, such that - unlike Eq. (fTT]) - the 
complexity of Pn{t) will increase with n. 

(iii) Under the same (infinite and high bias, respec- 
tively) assumptions we may adjust bias voltage and tem- 
perature such that we can distribute the left-associated 
Fermi-functions in an approximately equidistant man- 
ner between zero and one (such as e.g. fL{^d) ~ 1, 
hicd + U) = 2/3, fd^d+m = 1/3, and fLi^d+iU) = 
for N = 3), wc can analytically extract the current 
((7T,^(t))) = £ and the long-term scaling of the next 
higher cumulants (for 7^ 0) of P„(t) 



r (mif + 3N7]r + 

16iV?7 (r -f 77) 
T_ L ^ 3r^(2r + 377) 
64 7]N{T + rjf 



t. (12) 



These expressions demonstrate that the higher cumulants 
diverge for small 77 in the long time limit. In the limit 
of an infinite number of multi-stable channels {N — > cxo) 
however, this divergence is overshadowed by the expo- 
nentially large degeneracy of intermediate currents: If 
the control jump matrices in Eq. Q did not scale with 
A'^, the divergence of all higher than second cumulants 
would persist also in the limit TV — > 00. 

(iv) Without these assumptions, we can still nu- 
merically perform both the inverse Laplace transform 
of M{x,z) and a Fourier integral to obtain P„(i), 
which is typically evaluated using the saddle-point 
approximationi^ [compare also inset of Fig. (Hb)]. The 
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FIG. 3: (Color online) Probability distribution Pn(t) for the 
number of tunneled particles n at different times t (orange 
region). For increasing times [(a) — >■ (d)] the average of 
the distribution moves linearly towards larger n. It pro- 
ceeds a crossover from multi-modal to unimodal with a tran- 
sition time of r]~^ and becomes (nearly) a broad Gaussian 
for t ^ rj~^ (d). For i] — even the long-term distribution 
depends on the initial mixture (symbols) here chosen as ana- 
lytic continuation of p to 77 = 0, which reflects in the different 
peak weights. For comparison: Distributions Pii'\t) for in- 
dividual currents Ik (blue, bold green, bold red, and black 
curves in the background). The saddle-point approximation 
(dashed lines) only captures the long-term behavior t ^ r)~^ . 
Parameters: N = 3, = 1, F = 0.1, r) = 0.0001, U = 1, 
/3 = 1, V" = V* = 5, y 00. 



result in the multi-stable bias regime of interest is shown 
for different times in Fig. |31 Choosing p for r/ > as 
initial mixture, the statistics is multi-modal with iV + 1 
maxima for t < 77"^ [Figs.[31[a) and (b)]. In contrast, the 
statistics becomes uni-modal for t > rj~^ [Figs. [3jc) and 
(d)]. In the long-term limit, the distribution for 77 > 
evolves essentially into a broad Gaussian (d). This is 
a general property of systems with cumulants linearly 
evolving in time. In contrast, when 77 = 0, Pn{t) will 
depend on the initial configuration for all times. For 
example, when one initializes in one of the subspaces 
k = 0, . . . ,N, one will observe the distributions Pi'^'' (t) 
for the individual currents (curves in the background 
in Fig. 13]). Starting in a statistical mixture for 77 = 
yields a multi-modal distribution even in the long-term 
limit (symbols in Fig. [3]) and leads to a divergence of all 
higher-order cumulants. 



C. Experimental parameters 

For the observation of a multi-modal distribution (e.g. 
bi-modal in Ref J^) the measurement time must lie be- 
tween F~^ and 77"^ . The rates in Ref.— are of the order of 
F « IkHz and 77 « IHz, respectively, such that the time 
of measurement can be estimated between 1 ms and 1 
s. For a distance of a hundred nm between transport 
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and control system, the Coulomb interaction strength in 
GaAs can be estimated to C/ w 1.2meV. Provided the 
picture of a single transport level is still valid (i.e., for a 
significantly larger on-site Coulomb interaction energy), 
pronounced multi-stability should be observable around 
temperatures of T « 14K. Larger distances or screened 
Coulomb interactions would lead to lower temperatures. 

V. CONCLUSIONS AND OUTLOOK 

We have studied stochastic multi-stable transport in 
terms of an n-resolved Master equation with a simple 
and scalable block-tridiagonal Liouville super-operator 
([3]). Multi-stability can be revealed by the full counting 
statistics even when the first moments are insensitive: 
For measurement times smaller than the switching rate 
between the distinct transport channels the distributions 
are multi-modal when the initial state is a mixture of 
the multiple steady states for jy = (this is the case for 
p for finite rf). This enables direct access to the num- 
ber of decoupled non-degenerate subspaces. For longer 



times or degenerate subspaces this is not possible. How- 
ever, unusually large higher-order cumulants may point 
towards intrinsic multi-stability. In case of sufficiently 
low temperatures multi-stable distributions may become 
effectively bi-stable. 

However, if the initial mixture is strongly localized 
within one of the multiple subspaces (this would be a 
consequence of a projective measurement), a particle- 
number measurement would result in the associated cur- 
rent with high probability and all other currents with low 
probabilities. Consequently, a sequence of repeated mea- 
surements would yield the switching dynamics observed 
for single-charge counting detector a^^'^^i^^ . 

We finally remark that multi-stable behavior can also 
emerge due to the effect of coherences'^. 
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